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The "Color Glass Condensate" is an effective theory description for the 
small momentum fraction x degrees of freedom in a high energy hadron or 
nucleus, which can be understood in terms of strong classical gluon fields. 
We discuss the resulting picture of the initial conditions in a relativistic 
heavy ion collision. We describe recent work to show that the leading 
logarithms of the collision energy can be factorized into the renormalization 
group evolution of the small x wavcfunction. We then describe how this 
framework can be used to understand the long range rapidity correlations 
observed by the RHIC experiments. 

PACS numbers: 13.85.Hd,24.85.+p,25.75.-q 

1. The little bang of an ultrarelativistic heavy ion collision 

Quark gluon plasma is studied in the laboratory in collisions of heavy 
nuclei at ultrarelativistic energies, presently = 200^ GeV at RHIC in 
Brookhaven or in the near future 5500^ GeV at the LHC in CERN (where 
A is the atomic number of the nucleus). The collision process is a compli- 
cated one, starting from the formation and equilibration of the matter to 
its evolution in time and space and ending in the decoupling of the system 
into the hadrons that are observed in the detectors. 

The typical transverse momentum scales of the bulk of particles pro- 
duced is in the GeV range, much less than the collision energy. Thus 
the initial conditions depend on the small x ~ VTl\fs < 0.01 part of the 
nuclear wavefunction. Because of the In 1 jx enhancement of soft gluon 
bremsstrahlung this is a dense gluonic system. When the occupation num- 
bers of gluonic states in the wavefunction become large enough, of the order 
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of 1/as (meaning that the gluon field is of order 1/g), the nonhnear 
interaction part of the Yang-Mills Lagrangian becomes of the same order 
of magnitude as the free part. The relevant comparison is between the two 
terms in the covariant derivative = + igA^: the momentum scales 
= —id^ < gAf^ become nonlinear. In the small x wavef unction the rel- 
evant component is the transverse momentum, we are therefore led to the 
concept of a transverse momentum scale Qs, the saturation scale, below 
which the system is dominated by nonlinear interactions. When the colli- 
sion energy is high enough (x small enough) , Qs S> ^qc^^ and the coupling is 
weak: we are faced with a nonperturbative strongly interacting system with 
a weak coupling constant. On the other hand, the large occupation numbers 
mean that the system should behave as a classical field. This suggests a way 
of organizing calculations that differs from traditional perturbation theory. 
Instead of developing as a series of powers in gA^ we want to calculate the 
classical background field and loop corrections (which are suppressed 
by powers of g) to all orders in gA'^^ . The classical gluon field will then be 
radiated by the large x degrees of freedom, which we shall treat as effective 
classical color charges. This picture of the high energy wavefunction is re- 
ferred to as the Color Glass condensate (CGC, for reviews see e.g. [I]). The 
collision of two such systems leads, in the early stages l/\/s <^ r < ^/Qs, 
to classical field configurations known as the Glasma [2j. 

At early times (r <C Rj^, see Fig. [T] for the coordinate system) the bulk of 
the system cannot, by causality, be aware of its finite size in the transverse 
plane. It will therefore be in a longitudinally expanding, to a first approx- 
imation boost invariant (9^ = 0) state, a 1-dimensional Hubble expansion. 
Boost invariance can come in two flavors. As we shall argue in the follow- 
ing, the very early time glasma degrees of freedom are boost invariant at the 
level of field configurations. This means that the longitudinal momenta of 
particles redshift towards zero ~ l/r while px ~ constant and the system 
becomes very anisotropic in momentum space. This field level invariance 
is broken by quantum fluctuations suppressed by Og, which then eventually 
evolve into a more equilibrated fluid that is isotropic its local rest frame. 
What remains is a boost invariant proflle of particle flow, as in the Bjorken 
hydrodynamical picture. What concerns us in this paper is the very earliest 
glasma stage and the initial quantum fluctuations that serve as the seeds of 
isotropization. 

In the following we shall first discuss the leading order, classical field 
level, results for the structure of the glasma fields and gluon production. 
We shall then, in Sec. [3l describe some ingredients of the recent proof [3] 
that shows how the leading logarithmic divergences of the NLO corrections 
to these fields can be absorbed into the renormalization group (RG) evo- 
lution of the weight functionals describing the hard sources, the JIMWLK 
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Fig. 1. Spacetime structure of the CGC and glasma fields. It is convenient to use 
the coordinate system with proper time r = v2ar^z?+ and the spacetime rapidity 
1] = ^ \nx'^ /x~ . 

factorization theorem. In Sec. H] we shall then describe multigluon correla- 
tions in the same framework. 

2. Gluon production to leading order and the glasma 

The CGC framework is based on a separation of scales between small x 
and large x degrees of freedom, which are treated as a classical field and an 
effective color charge density. In practice the classical field is obtained from 
the equation of motion 

[D^,F'^n = J"- (1) 

The current in the case of a nucleus-nucleus consists of two infinitely Lorentz- 
contracted (this picture will be discussed more below) nuclei on the light 
cone [1]: 

J'' = <5'^+P(i)(^±)5(x-) (2) 

The large x degrees of freedom have now been reduced to a classical effective 
color charge density p{x±), which is a static (hence the "glass") stochastic 
variable. Its values are drawn from a probability distribution Wy[p{x±)] 
which depends on the cutoff rapidity y = In 1/x separating large and small 
X. To a first approximation we can take e.g. the Gaussian distribution of 
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Fig. 2. Left: Gluon spectrum from the leading order classical field computation. 
Right: the components of the glasma field, the initial condition is a longitudinal 
electric and magnetic field, the transverse components develop in a time ~ 1/Qs- 



color charges that defines the MV [5] model 



W[p{x±)] = AAexp 



(3) 



The probability distribution Wy[p{x±)] is analogous to a parton distribution 
function in the DGLAP formalism; it is a nonperturbative input that we 
are not able to compute from first principles, but one can derive evolution 
equation for its y-dependence. This equation is known by the acronym 
JIMWLK. 

For a fixed configuration of the color sources p the calculation of the 
Glasma fields proceeds as follows P] . The solution of the Yang-Mills equa- 
tions in the regions of spacetime > 0, < that are causally connected 
to only one of the nuclei (areas (1) and (2) in Fig. [1]) is an analytically known 
pure gauge field. It gives the initial condition for the numerical solution in 
the forward light cone (3). Working in the temporal gauge = these 
initial conditions are 



^^=0 = ^il)+^(2) (4) 
^{2)J' 



where 2) pure gauge fields that are the solutions of the one- 

nucleus problem 

^ii,2)='-U(i,2){x^)d,ul,^{x^). (6) 
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These pure gauge fields are gauge transforms of the vacuum with the Wilson 
lines computed from the color charge density 

t/(i)(:rx,x-)=Pexp|-i5 y"dy-^^^^^|, (7) 

with the Wilson line ?7(2) given by the analogous formula in terms of the 
other color charge density. 

These Wilson lines are in fact the most natural variables to describe the 
soft gluonic field degrees of freedom of the nucleus; they correspond to the 
eikonal scattering amplitude of a color charge off the strong color fields. For 
example the dipole cross that determines the structure function measured 
in deep inelastic scattering is a correlator of these same Wilson lines. The 
upper limit of the ^"-integral in Eq. ([7]) must be thought of as x~ ~ e^; 
when the cutoff rapidity y becomes larger (x smaller), smaller momentum 
p'^ gluons are considered as part of the source, which consequently extends 
further in the conjugate variable x~ . Thus each infinitesimal step in the 
renormalization group evolution towards smaller x corresponds to adding a 
layer in x~ to the color source, or equivalently to multiplying the Wilson 
line by an SU(3) matrix that is infinitesimally close to identity. 

From the point of view of the classical glasma fields in Eq. ^ the Wilson 
lines are independent of the longitudinal coordinate: the longitudinal struc- 
ture appears only indirectly in the properties of the probability distribution 
Wy[C/]. The classical fields represent degrees of freedom with a smaller p"*" 
than the ones integrated out to the Wilson lines and are not able to resolve 
their structure which is shorter range in x~ . This is the sense in which the 
(5-functions in the currents of Eq. ([2|) must be understood. 

The numerical method for solving the Yang-Mills equations in the for- 
ward light cone was developed in Ref. [6j and the actual computations re- 
ported in Ref. [7] The equations of motion are most conveniently solved 
in the Hamiltonian formalism. Due to the boost invariance of the initial 
conditions in the high energy limit the Yang-Mills equations can be dimen- 
sionally reduced to a 2-1-1 dimensional gauge theory with the ry-component 
of the gauge field becoming an adjoint scalar field. With the assumption of 
boost invariance one is explicitly neglecting the longitudinal momenta of the 
gluons. In the Hamiltonian formalism one obtains directly the (transverse) 
energy. By decomposing the fields in Fourier modes one can also define 
a gluon multiplicity corresponding to the classical gauge fields; the result- 
ing gluon spectrum is shown in Fig. [2j The color fields of the two nuclei 
are transverse electric and magnetic fields on the light cone. The glasma 
fields left over in the region between the two nuclei after the collision at 
times I < T < 1/Qs are, however, longitudinal along the beam axis [2] (see 
Fig. ED. 
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Fig. 3. Factorization of LLog corrections to gluon production: the phase space in- 
tegral over Ay diverges and is cut off at the separation scales 2/1,2- The dependence 
of the color charge density distributions Wy-^ ^ on the cutoff cancels the leading 
logarithmic part of the dependence on y\^i- 



3. Factorization 

To understand the context of the high energy factorization theorem 
proven in Ref. [3j it is perhaps useful to look first at the weak field limit of 
the CGC, where particle production can be computed using /cy-factorization 
([8|, see e.g. [9] for an application to heavy ion collisions). The leading order 
multiplicity is 

dA^ 11/" d2fc± ^ , 

To obtain the real part of the leading log correction to this result one must 
take the corresponding expression for double inclusive gluon production 

dAf 11/" d2fc_L , , , s . N 

and integrate it over the phase space of the second gluon {q±,yg)- Note 
that at leading log accuracy we have here taken the multi-Regge kinematical 
limit, assuming that the two produced gluons are far apart in rapidity (see 
e.g. [TO]). The integral over yq diverges linearly (this is the general behavior 
of the gg — > gg scattering amplitude in the high energy limit t fixed, s ~ 
—u — > oo). This divergence is compensated (to the appropriate order in Os) 
by the real part of the BFKL evolution equation for ipy(k±i). 

In the fully nonlinear case of AA collisions the fcj-factorization is broken 
(see e.g. [6l[TT]), and one must solve the equations of motion to all orders in 
the strong classical field. The analogue of the unintegrated parton distribu- 
tion (py{k±) is the color charge density distribution VFj^[p]. These are similar 
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in the sense that they are not (complex) wavefunctions but (at least loosely 
speaking) real probability distributions. Factorization can be understood 
as a statement that one has found a convenient set of degrees of freedom in 
which one can compute physical observable from only the diagonal elements 
of the density matrix of the incoming nuclei. The difference is that when 
in the dilute case these degrees of freedom are numbers of gluons with a 
given momentum, in the nonlinear case the appropriate variable is the color 
charge density and the relevant evolution equation is JIMWLK, not BFKL. 
The kinematical situation, however, remains the same. To produce a gluon 
at a very large rapidity (or a contribution in the loop integral of the virtual 
contribution with a large k~^) one must get a large +-momentum from the 
right-moving source. Thus one is probing the source at a large , i.e. small 
distances in x~, and the result must involve Wy[p] at a larger rapidity (see 



The underlying physical reason for factorization is that this fluctuation 
with a large requires such a long interval in to radiated that it must 
be produced well before and independently of the interaction with the other 
(left moving and thus localized in x^) source. The concrete task is then to 
show that when one computes the NLO corrections to a given observable in 
the Glasma, all the leading logarithmic divergences can be absorbed into the 
RG evolution of the sources with the same Hamiltonian that was derived by 
considering only the DIS process. This is the proof ^3| ll2t [T3] of factorization 
that we will briefly describe in the following. 

Consider the single inclusive gluon multiplicity which is a sum of prob- 
abilities to produce n + 1 particles, with the phase space of the additional 
n must be integrated out 



Because we have a theory with external color sources of order p ~ 1/g, 
all insertions of the sources appear at the same order in g [Ti]. A calculation 
using the Schwinger-Keldysh formalism leads to the following results: At 
LO, the multiplicity is obtained from the retarded solution of classical field 
equations (here (...) includes the appropriate normalization and projection 
to physical polarizations) 



The NLO contribution includes the one loop correction to the classical field 
and the H — component of the Schwinger-Keldysh (SK) propagator in the 
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background field 



diV. 



(12) 



Now consider a small fluctuation a^{x) of the gluon field around the classical 
value. The H — (SK index) component of the propagator is bilinear in these 
small fluctuations satisfying retarded boundary conditions. Also the virtual 
term [5 satisfies an equation of motion with a retarded boundary condition 
and a source term involving a loop in the classical background field, see 
Fig. [3] for a pictorial representation of this structure. One can express the 
propagation of such a small fluctuation a^{x) above the past light cone S 
as a functional derivative of the LO classical field A^{x) with respect 
to its initial condition on S: a^{x) = J^g^ a{'il) ■ TuA^{x). This leads after 
some rearrangements to the expression for the NLO contribution to the 
multiplicity as a functional derivative operator acting on the leading order 
result: 



dN 






d^p 


NLO 


2L 



dN 
d^p 



■ (13) 



This expression involves the part of the two point function below the light 
cone S: 

G'^i^,^) ^ / ./'L at^{u)a\^{v). (14) 



(27r)32Sfc 
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Here the small fluctuation field a^{x) is the solution of the linearized equa- 
tion of motion in the classical field background with an initial condition 
given by a plane wave lim^o ^ _ a'^f,{x) = e^{k)e^'^^'^ . 

The leading logarithmic contribution comes from the longitudinal com- 
ponent of the integral over fc, the momentum of the initial plane wave per- 
turbation (and the corresponding momentum in the one loop source term 
for the equation of motion satisfied by (5) . This LLog part of the functional 
derivative (fT3]) operator turns out to be precisely equivalent to the sum 
of the JIMWLK Hamiltonians describing the RG evolution of the source 
distributions The JIMWLK Hamiltonian 

is most naturally expressed in terms of Lie derivatives Da{xj_) operating 
on the Wilson lines introduced in Eq. ([7]). in terms of which the kernel in 
Eq. (dl]) is 



VT y {x± - u±y{yj_ - u±y i 



ab 



-Uix^)U\u^)-U{u^)U\y^) + l\ . (16) 

The fact that no other terms with the same logarithmic divergences appear 
is the proof of factorization; this is the central result of Ref. [3]. 



4. Multigluon production 
4.1. Short range in rapidity 

Let us then consider the probability distribution of the number of gluons 
produced in a small rapidity interval. It was shown in Ref. yL2j that a similar 
factorization theorem holds for the leading logarithmic corrections to this 
probability distribution in the sense that we will briefly review here. It is 
convenient to deflne a generating functional 

00 

n=0 

The Taylor coefficients of J- around z = 1 correspond to the moments of 
the probability distribution; integrated over the momenta of the produced 
gluons they are 



d"iV„ 



d3pi • • • d3p^ 



(17) 



{N) {N{N-l)) ... {N{N -1)---{N -n + l)) . (18) 
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The result of Ref. [12] is that when these moments are calculated to NLO 
accuracy, the leading logarithms can be resummed into the JIMWLK evo- 
lution of the sources completely analogously to the single inclusive gluon 
distribution. The resulting probability distribution can be written as: 

d"P„ f ... r .... r n 1 dN dN 



d^pi • • • d3p„ 

P1.P2 

Note that the Poissonian-looking form of the result is to some extent an arti- 
fact of our choosing to develop and truncate precisely the moments Eq. (jl7p 
that are simply (A^)" for a Poissonian distribution. Since in our power 
counting N ~ l/ctg, any contributions that would make the distribution 
Eq. (jl9p deviate from the functional form are of higher order in the weak 
coupling expansion of the moments (|17p and are neglected in our calculation 
unless they are enhanced by large logarithms of x. Nevertheless it should be 
emphasized that in spite of appearances of Eq. (|19p the probability distri- 
bution is in fact not Poissonian. To understand the nontrivial nature of this 
result it must be remembered that the individual factors of in Eq. (fT9]) 
are all functionals of the same color charge densities pi^2] thus the averaging 
over the p's induces a correlation between them. These correlations are pre- 
cisely the leading logarithmic modifications to the probability distribution; 
they have been resummed into the distributions Wy] the functional form of 
the multigluon correlation function under the functional integral in Eq. ()19p 
is the same as at leading order. This is the result of the proof in Ref. [12]. 



4-2. Long range in rapidity 

Let us now relax the restriction that the gluons should be observed 
only in a small rapidity interval and allow for arbitrary separations in ra- 
pidity [15j. There is now phase space available to radiate gluons (and for 
the corresponding virtual contributions) between the measured gluons, and 
including this radiation can introduce additional large logarithms of the 
energy. To develop a physical picture of this situation it is perhaps use- 
ful to take a step back and consider a more general picture of JIMWLK 
evolution in terms of its Langevin formulation derived in Ref. [TB]. The 
original derivation is presented purely as an alternative formulation to gen- 
erate the single Wilson line probability distribution that solves the JIMWLK 
equation. The JIMWLK equation as it is usually written, as an equation 
satisfied by the probability distribution of Wilson lines at a single rapidity 
y, does not formally give information about correlations between differ- 
ent rapidities. Going back to the derivation one sees, however, that the 
rapidity correlations are also encoded in the formalism. This is most trans- 
parent in the Langevin formulation, where one can identify each trajectory 
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in the Langevin equation with one high energy colhsion event. In this sense 
the Langevin formulation contains more physical information than just the 
JIMWLK equation for the probability distribution at a single rapidity; it 
also gives the combined probability distribution for Wilson lines at different 
rapidities 

Wy,...yAU,{x^),...,Un{x^)]. (20) 

Knowing the general (multiple rapidity) probability distribution will enable 
us to compute the correlations between Wilson lines, and consequently of 
physical observables such as the multiplicities, at different rapidities. In 
the following we will give a more precise formulation of this statement and 
show that it is consistent with our previous result concerning multigluon 
production. 

In the Langevin formulation the distribution W of the Wilson lines can 
be obtained by evolving in rapidity the elements of an ensemble of Wilson 
lines according to 

U{y + dy, x^) = U{y, ;r^)e-*"(^'^^)'^^ (21) 

where the change in a small step dy in rapidity is given by a deterministic 
term and a stochastic term, 

a^{y,x^) = a-{xuy)+ [ e"^(^x, zx)C'(zx, y). (22) 

J Zj_ 

Here is a Gaussian random variable defined by {Ci{x±,y)Cj{y±Ty')) = 
5°''^5ij6'^{x± — y±)S{y — y') and the square root of the JIMWLK kernel is 

e'^^^L.z^) ^ }^ -'^~^\^ {l-U\x^)U{z^)r. (23) 

This stochastic formulation is the method used in numerical studies of the 
JIMWLK equation [17]. 

Knowing that the correlation follows from a Langevin equation imposes 
an additional structure (of a Markovian process) on the probability distri- 
bution: 

Wy,,y, [U^, U'^] = Gy^^y^ [U^, U^] Wy^ [C/^] , (24) 

where the JIMWLK Green's function G is determined by the initial condi- 
tion 

hm Gy^^y^ [U^\ U'^] = 6 (U^ix^) - W^ix^)) (25) 

Vp^Vq 

and the requirement that it must satisfy the JIMWLK equation 

dypGyp-yq P", U^] = n {U^{x^)) Gy^^y^ [U^, . (26) 
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This JIMWLK Green's function contains all the information, at the leading 
log level, of long range rapidity correlations in gluon production. This 
structure follows from the computation of the leading log part of 1-loop 
corrections to a wide class of observables that can be expressed in terms of 
correlators of the gluon fields at r = (or equivalently Wilson lines) 



{0)^^^^ = j [DUi{y,x^)] [DU2{y,x^)]W[U^{y,x^)]W[U2{y,x^)]0,^. 

(27) 

Here we have introduced a continuous rapidity notation W\U{y, a3_|_)]for the 
probability distribution of Wilson lines (j20p . This should be understood as 
a probability distribution for the trajectories that the Wilson line U{x^) 
takes on the group manifold along its evolution forward in y following the 
Langevin equation. We can formally return from the distribution of trajec- 
tories to a distribution of Wilson lines at one individual rapidity as 

Wy[U{x^)] = j [DU{y, x^)] W [U{y, x^)] 5 [U{x^) - U{y, x^)] . (28) 

Equation (I27p is the central result of Ref. [15] , showing that all the lead- 
ing logarithms of rapidity (either the rapidity intervals between the nuclei 
and the tagged gluons, or between the various produced gluons) can be ab- 
sorbed into the probability distributions W for the trajectories of Wilson 
lines of the two projectiles. However, the crucial point to keep in mind is 
that it involves an average over y-dependent "trajectories" of Wilson lines, 
rather than an average over Wilson lines at a given fixed rapidity. 

The calculation of multigluon correlation is in fact simplified in the 
strong field limit, where the leading contribution to particle production 
corresponds to the classical field and the correlations are encoded in the 
evolution of the sources [I8j . In the "pA" case where one of the sources is 
assumed to be dilute, the situation becomes much more complicated, be- 
cause the disconnected classical contributions are not the only dominant 
ones any more. This structure is illustrated in Fig. O 



4-3. Application to multigluon correlations 

Let us now specialize Eq. (j27|) to the case of the single and double 
inclusive gluon spectra. The single inclusive gluon spectrum dNi/d^p at 
LO depends only on Wilson lines C^i,2(yp5 ic±) at the rapidity yp of the 
produced gluon. One then obtains the known result for the single inclusive 
gluon spectrum as 



dA^i 



d2p±dy 



LLog 



[DUl] [DU2]Wy^ [Ul] Wy^ [U2 



dNi[Ui,U2 



d2p_Ldy 



(29) 
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Fig. 5. Relative importance of connected and disconnected diagrams to the two 
gluon correlation function. One of the color charge densities is considered large, 
Pi ^ whereas the other is allowed to vary between the "AA" case P2 ~ l/ff 
and the "pA" one p2 ~ 5. The order of the disconnected diagram, on top, is 
g'^pfpi, whereas the interference diagram in the middle is g^PiP2 and the connected 
one, lowest, is g'^p\p2- In the "AA" case the disconnected diagram dominates, for 
the "pA" case all three are equally important. In the dilute "pp" limit only the 
connected diagram matters and both gluons are produced from the same BFKL 
ladder. 



For the resummed inclusive two-gluon spectrum, we must recall that at 
LO it is simply the product of two single gluon spectra, each of which 
depends on Wilson lines at the rapidity of the corresponding gluon. It is 
then straightforward to proceed as in the case of the single gluon spectrum 
in order to obtain : 



dN2 



d2p_Ldypd2g_Ldyq 



LLog 



[DUl] [DUl] [DUl] [DUl] X 



(30) 



where the double probability distribution Wy^^^y^ [?7f , U^] is given by Eqs. ([2 
and (ESI). 
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We have now assembled all the ingredients needed to compute rapid- 
ity correlations in the Glasma and address features such as the elongated 
"ridge" structure in the two particle correlation observed at RHIC |19j . 
There have already been several boost invariant classical field calculations |20j 
of this effect and the azimuthal structure, but the inclusion of quantum evo- 
lution is needed to understand the rapidity dependence. As a first approx- 
imation one should be able to formulate the equivalent of the mean field 
approximation leading to the BK equation for the JIMWLK propagator. 
Numerical studies of the JIMWLK equation would then be needed to study 
the validity of this approximation for rapidity correlations; in the structure 
of the single rapidity probability distribution the violations from the mean 
field limit have been observed to be small [21j . 
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